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SECTION 1.0 
INTRODUCTION 


In 1973, the Radio Astronomy Explorer-B (RAE-B) 
satellite will be placed in a 1100 km-altitude circular 
orbit around the Moon to make radio astronomy measurements. 

The purpose of the simulation study described in this 
report is to determine the 136 MHz and 400 MHz noise temper- 
ature of the ground network antennas which will track the 
RAE-B satellite during data transmission periods. Since the 
.noise temperature of the antenna effectively sets the Signal- 
to-Noise Ratio (SKR) of the received signal, a knowledge of 
SNR will be helpful in locating the optimum time windows for 
data transmission during low-noise periods. 

Antenna-noise temperatures at 136 MHz and 400 MHz will 
be predicted for selected earth-based ground stations which 
vj^ill support RAE-B. Telemetry data acquisition will be at 
400 MHz; tracking support at 136 MHz \tfill be provided by the 
Goddard Range and Range Rate (RARR) stations . 

The antenna-noise temperature predictions will include 
t]ie effects of galactic-brightness temperature, the sun, and 
the brightest radio^ stars. Predictions will cover the ten- 
month perj od from March 1, 1973 to December 31, 1973. The 
RAH-B mission A\-ill be especially susceptible to SNR degrad- 
ation during the two eclipses of the Sun occuri'ing in this 
period . 


The body of this report is divided into three parts: 

o Development of Equations 

« Computer Progr.nm Description 

c V.'orst-Casc Analysis Results 
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The worst-case analysis is a "first-look" approx- 
imation to determine preliminary estimates of antenna-noise 
temperature profiles. A more refined computer program is 
being written that predicts antenna-noise temperature more 
accurately by utilizing complete models of the antenna patterns 
and station location geometry. The refined model is based 
upon the Satellite Data Quality Program previously developed 
for GSFC/NASA under contract Number NAS 5-11736 PCN-523-W- 
70446. 



SECTION 2.0 
PROGRAM DESCRIPTION 


2.1 DEVELOPMENT OF EQUATIONS 

Four sources o£ antenna-noise temperature are con- 
sidered in this study: 

• ■ Sky-Brightness Temperature 

• Sun 

• Radio Stars 

• Antenna Back Lobe Noise Temperature 

o Total Antenna-Noise Temperature 

The formulation utilized within this program has 
been previously utilized in the -Data Quality Prediction 
Program [1] developed by Wolf Research ^ Development Corp. 
for NASA/GSFC under contract NAS 5-11736 DCN 523^W-70446. 
The equations presented in the follov/ing are taken from 
References 1 and 2 . 

a) Sky-Brightness Temperature 

Kraus [3] develops the following formulation for 
sky briglitness temperature 


e=9O°-0o <|>=2Tr 

/ T (6,<j») G (0,<|)) sin 0303$ 

n *4l 


"^SKY* 


0»=9O“-0o 

fo Jo 


G (0,<^) sin 0303(J) 
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wlicrc 


0 

0 ■ = elevation angle between antenna's boresight 

axis and the horizon, degrees 

is the noise temperature distribution of the 
galaxy (excluding the sun and predominant 
radio stars) obtained from References 1 and 2. 

G(6,<)) is the lossless antenna gain distribution. 

9 and (}> are the orientation angles defining the 

position of a radial surface element within 
the celestial hemisphere. 

Figure 1 shows the physical relationshio of the 
variables given in the above equations. The double integral 
is computed by rectangular integration assuming that the 
antenna boresight is directed at the center of ‘ the Moon's 
optical disk. The antenna patterns and brightness tempera- 
tures are accessed from magnetic tape storage in the computer 
program. 
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Rcloficn. of Antcnnc f'ollcrri to Ocicttiol 


Sphorc 


Figure 1. 


This program required the use of an accurate radio-sky 
map which covers the celestial sphere completely, for both 
136 Mllz and 400 MHz. Since detailed sky-brightness temperature 
contours were not available at these frequencies, it became 
necessary to scale either existing radio maps in temperature, 
or to generate a composite map from various smaller maps. The 
136 MHz and 400 MHz sky-brightness temperature maps, appearing 
in Appendix A- 3, were prepared in this manner. 

The 136 MHz radio map was scaled from data at 150 MHz, 
publislied in 1971 by Landecker and Wielebinski (Reference 4). 
The following relationship was used for scaling: 



degs. K. 


Reference 5 was utilized for the conversion of galactic 
coordinates, employed by Landecker and Wielebinski, into the 
necessary equatorial coordinates required for this program. 

The 400 MHz radio map, appearing in Appendix A- 5, is 
a composite map formulated from tlie sectional maps published 
in 1962 by Pauliny-Toth and Shakeshaft (Reference 6), and 
in 1956 by Droge and Priester (Reference 7) . 


b) Sun 

The contribution of the sun to the antenna-noise 
temperature is given by Berkowitz [8] as: 



T, , assuming G >> 0 
o' “a s 
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v:herc 


0 = Angular radio dinweter of sun*s apparent temper 

aturc model, degrees; assume * 0.66 at 136 
I MHz and 400 MHz. 

0 = Half-power beam width (HPBW) of symmetrical 

3 . 

antenna main lobe. 

T|^ = 8 X 10^ K for quiet sun ideal model at 136 MHz, 

and. 6 X 10^ K at 400 MHz. 


c) Radio Stars 


The following equation is used by Taylor [1,9] to 
compute antenna noise power rise due to a point- source radio 
star: 



n=l 


G • G(0) 

p LI , 

4tt 


Do 


for a single polarization. 


where 

M - is the number of radio stars 

p - is the observed radio star noise flux density , 
wm’^Hz’^, constant over band\Hdth Af. 


7 



X - wave length of transmission 

Gp - ]u\'!h antenna power gain, above isotropic 

G(G)- antenna gain attenuation at angle 0 off-boresight 
i.e., G(0)=1 for 0 *0. 

Note that the antenna-noise temperature is 
N* 

T* * degs . K 

kA£ 

where 

Njs = total noise power. due to all radio stars 
within the antenna's radiation pattern 

- 23 

k = Boltzmann's constant, 1.38x10 J/K 
A£ «. noise bandwidth of receiver, Hz 


d) Antenna Back Lobe Temperature • 

The black-body radiation of the Earth contributes to 
the overall antenna-noise temperature by means of the back 
lobe of the ground antenna, 

*. 

Based on Blake's data (Reference 10), the effect of 
antenna back lobe temperature, is approximated by 

adding a constant to the equation for antenna-noise temper- 
ature ns follows: 
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136 Mllz 


"^BACK 




400 Mllz 


^^BACK 


35°K 


c) Total Antenna-Noise Temperature 

The total antenna-noise temperature, T^q.j, is computed 
by siiTTiiP-ing each of the four contributions • 


Ttot 


■^SKY "^SUN '^BACK 


Antonna-noise temperature will be maximum at New Moon, 
once each month. A higher peak is reached once each year 
CDecember) when the Galactic Nucleus is eclipsed by the Sun, 
during Kev/ Moon (see Figure 2) . 




tL< 


10 


I 


gure 2. Celestial Source Spatial 
Arrangement in Month of December 



2.2 


PROGRAM LOGIC PLOK 


The follov.’ing logic steps determine the antenna-noise 
tcriporaturc for cacli ground antenna Khich will track the 
RAR-B; 

• Each station is tested for a visible moon. 

6 Antenna-noise temperature is computed for 

each station having a visible moon. 

o Antenna-noise temperature is recomputed at 
periodic time intervals from Moonrise to 
Moonset. 

The equations and logic for each major portion of the 
program will now be described and follov/ed by a Program flow 
chart. 
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2 . 2.1 


Visihnity Tost (Block 1 to 4 In Snction 2.2.5 ) 


Tor each .statioi, the zenith vector is computed 


.TP ; 



cos <|> cos X 
cos <j) sin X 
sin ^ 


where 


(J> is the geodetic latitude 
X is the geodetic longitude. 


Klien the dot product of the moon's unit vector, T... and T^, is 
negative the station is located below' the horizon and the 
antennn-noise temperature is not calculated. Tj^j is computed a 

V “= ^m" 
cos sin 

sin (j) 


where 



declination of moon at the time antenna-noise 
temperature is computed 


and 


Xjjj is the longitude of the moon. 



U'hcn tlie dot product of Xp and Tj^j is not negative the station 
antonn.'i-noisc temperature is calculated. This visibility test 
is repeated for every station vithin the same time period. 


2.2.2 Antenna -Noise Temperature Calculation (Block 5 to 16 

In .Ejection 2.2.5) ” 

a. Contribution of Sun to Antenna-Noise 
Temperature (Block 5 to 10) 

In order to determine if the sun's radio diameter 
falls v.’ithin the antenna main lobe, the sun's position unit 
vector (Tg) is computed in the earth- fixed coordinate system 
as 



cos <J)g cos Xg 
cos <{> sin A 
sin 


where 


(b is the declination of sun at a given time and 
s 

Xg is the longitude of the sun. ^ 

If the dot product of the moon's unit vector, Tj^j, and the 
sun's unit vector, Tg, is less than the half power beam width 
of the antenna's main lobe, then the antenna-noise temperature 
due to the sun is calculated as shown in Section 2.1. Other- 
wise, the effect of the sun is assumed negligible. 
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b; Contribution of Radio Stars to Antenna-Noise 
Tei/ju raturc (block 11 to 14) 

Tlie position unit vector of each radio star is computed 
in tlic c.'irth-f ixod coordinate system as 



COS 

*^SK 

CCS 

^SR 

T 

•^STAR 

cos 

^SR 

sin 

^SR 


sin 

^SR 




.where 


is the declination of radio star at a given 
time 



the longitude of a given radio star. 


When the arc cos of the dot product of the moon's unit vector, 
Tj^j, and radio star's usiit vector, is less than one-half 

the beam width of antenna main lobe, the computation is per- 
formed as shown in Section 2.1. Otherwise, the effect of a 
radio star is considered negligible. For each radio star, a 
computation is performed. The antenna-noise temperature due 
to all radio stars is the sum of the temperature for each 
radio star. 


c. Contribution of Sky Temperature To Antenna-Noise 
Temperature (Block 15) 

Sky-noise temperatures are stored as discrete values 
at regular intervals. The total noise power is a function 
of the noise-temperature distribution of the sky, in convolution 
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with the antenna j^ain distribution over the celestial sphere 
(sho’..a Scct.icm 2.]). The nuiucrical integration required 
for perforrieJ utilizing the shy map. 

d. Contribution of Antenna Back Lobe to Antenna- 
Noise Temperature (Block 16) 

Back lobe noise temperature is assumed constant at: 
Tb^CK “ ^ “^BACK " ^ 

e. Total Antenna-Noise Temperature (Block 16) 

-The tracking station's total antenna-noise temperature 
is the sum of contributions from the sun, radio stars, sky 
background temperature, and Antenna Back Lobe temperature. 


2.2.3 Terminal Test 


The above operations (Sections 2.2.1 and 2.2.2); are 
performed for each station antenna, for the same time interval. 
After all the stations have been tested for a visible moon, 
the total antenna-no,ise temperature is computed for the initial 
antenna-pointing angle; the computation is repeated at an 
incremental pointing angele, until Moonset. 


2.2.4 Output (Blocks 19 and 20, Section 2.2.5 ) 

Plots of antenna -noise temperature versus time, for 
each station, are output using the V/OLF Plot Package 
(Reference 11) . 
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compute the unit vector, of sun 


-dot product of sun's and 
moon's unit vectors 


--compute radio star's 
unit vector 


--dot product of moon's and 
radio star's unit vector 















SKCTION 3.0 

SlMPUFIl-P COMPUTATION- 01- ANTPNNA-NOISE TEMPERATURE 


3.1 APPROXIMATIONS AND ASSUMPTIONS 

•> 

The Simplified Computation program uses the logic and 
equations in Section 2 to predict the worst-case antenna-noise 
temperature for the ground station antennas. This program 
utilizes a number of approximations to obtain quick- look 
worst-case results. This program is based upon the following: 

1. Time span for observation is for the last ten 
months in 1973 (March 1 to December 31) . The 
time interval for computation is variable between 
10 minutes to 60 minutes, the time span per run 
being variable. 

2. Ephemeris tape (prepared by the Jet Propulsion 
Laboratory (JPL)) provides x,y,z coordinates 
(referenced to equinox) of Sun and Moon at any 
time. The station unit vector is computed from 
right ascension and declination, instead of 
geodetic latitude and longitude. 

3. Antenna radiation pattern is approximated by 

an ideal, rectangular, gain distribution corre- 
sponding to the half power beamwidth of the main 
lobe. Effects of antenna side lobes are consid- 
ered negligible with the exception that if the 
first side lobe impinges upon the Sun or the 
Galactic Nucleus, a uniformly- illuminated side- 
lobe of the appropriate gain is also employed 
in the calculation. 
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4. Center of the lunar disk is assumed aligned vjith 
tl ':0 .'intenna's borcsight axis. 

5, Program ignores antenna back lobe temperature. 


3.2 SIMPLiriKD COMPUTATION RESULTS 

3.2.1 Voorst Case Temperature Profiles for 1975 

The raaximuin antenna-nbise temperature is computed at 
136 MHz and 400 MHz for each station on every day between 
March 1 , 1975 and December 31 , 1973. This date is shov;n in 
Figures 3 through 15 which show the envelopes of the peak daily 
temperatures. The measurement interval in each case is 60 
miniites. Figures 3 through 8 plot the antenna-noise temp- 
erature at 136 MHz for tlie VHF Range and Range Rate 
antenna (UAl<r.) for all tracking stations supporting this 
antenna configuration. Figures 9 through 11 give the temp- 
erature profiles for the 85- ft. (26-m) diameter antenna 
at 400 MHz, and Figures 10 through 13 show corresponding 
data for the 40-ft. (12-m) diameter antenna at 400 MHz. 

The results show that a "cool sky" antenna-noise 
temperature for the 136 MHz antenna is approximately 500'K, 
while tlie corresponding value for the 400 MHz antennas is 
in the vicinity of 2S‘’K. 

Antenna-noisc temperature "hot spots" occur roughly 
at 29 j day intervals, corresponding to a New Moon (see Fig- 
ure 3), iv’hcn the Sun enters the antenna main lobe. The 
"Astronomical Phenomena For The Year 1975", issued by the 
Nautical .Almanac Office, United States Naval Observatory, 

'la.shi ngton , D.C., gi ves the Universal. Time (UT) for a New 
Moon (covering March 1 to December 31, 1973) as: 
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Kcw 
j'loon 
Hr! 1 '? 

March 5 

April 3 

May 2 

June 1 

June 30 

July 29 

August 28 

September 26 

October 26 

November 24 

December 24 


UT(hr. ,mln 0 

11^5’^ 

20 ^ 55 ”* 

04 ^ 34 ”* 

llhjgin 

18^ 59*" 
03^25"^ 
13^54'^ 
03 ^ 17 "' 
19^ 5 s’" 
15^07"* 


>- 
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MAR 1, 1973 TO DEC 31, 1973. TRACKING STATION; 



FTHirPR 
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KIAR 1 ..1973 TO DBC 31, 1373, TRACKING STATION; 





29 


I 
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973. .TRACKING STATION: 
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M\R 1 , 1973 TO nnc 31, 1973. TRACKING STATION 
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3.2.2 Fine-Grniin A nal ysis of Hiuli Temp erature Rcpions 

* Figures 16 thorugh 21 show fine-grain detail for 
typical high-tcinpcraturo regions. The zero temperature poitions 
of these curves represent periods when the moon is not visible. 
The peaks temperature recorded for each d.iy is utilized as 
data for the cunmulative temperature envelopes given in Figures 

i through 15. 
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RAE-B TRACKING ANTENNA TB*PERAT\JR£ 
TRACKING STATION » ORCRAL ANTEf4'(A TYPE * S5F 
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SECTION 4.0 
CONCLUDING REMARKS 


Tho I'osults of the simnlified simulation are expected 
to be accurate within tv;enty percent for tho stations em- 
ploying the VUF RARR or 40 foot dish antenna. However, the 
85 foot dish results may not be at the same .level of accuracy 
because the maximum error introduced by the geocentric sim- 
plification is almost the same magnitude as the half beam 
width of the 85 foot dish main lobe (1.25®). This simplifi- 
cation also causes a timing or phasing error in the temp- 
erature profile of betv^een one to two hours for all stations. 


These preliminary results will be verified and re- 
fined by the utilization of a more sophisticated program, 
now under development, which does not employ the approxi- 
mations of the simplified program. The temperature profiles 
obtained from this program will be published in a subsequent 
report. 

Appendices A and B have been inserted to further 
define the program. 


4.1 ANGULAR ERROR INTRODUCED BY GEOCENTRIC APPROXIMATION 

Consider the following geometry of the moon, earth's 
center and points of moon rise and moon set on the earth's 
surface 
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Angle ACD 

Angle ABD 

EC 

AC 

* 'I rBC-> 
tan (^) 

Thei'cf ore , 
approxinat 



A - moon's position 
B - moon rise 
C - geocenter 
D - moon set 


■= ACB = 90® since AC is the bisector of an 

isosceles triangle 

= ADB - 90® - tan“^ BC 

SC 

" Average Earth's radius 6378.155 km 

= Average distance from geccenter 
to’ selcnocentcr, 384000 km 

= 0.0166 = 0.95® 


the angular error introduced by the geocentric 
ion is a maximum of 0.95® 
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AP PUK’D IX A 
PROGRAM INPUT DATA 
FOR SIMPLiFIUD PROGRAM 


1. Time span for observation and interval of computation. 
For the last ten months in 1973 (from March 1, 1973 

to December 31, 1975). The interval of computation 
is variable betKccn 10 minutes to an hour. The span 
per run is variable, 

2. Name and location of the tracking stations are given 
in Section A-1. 

3. Antennas at 136 Mllz and 400 Mllz arc specified in 
terms of: 

a. Main Lobe half pov.'er beam width 

b. Second lobe total beam width and gain factor. 

The tracking antenna patterns are displayed in 
Section A. 2 . 

4 . Sky Map, Both 136 MUz and 400 MUz brightness tempera- 
ture maps are on tape which are accessed by the pro 
gram. Contour plots of these data are given in Section 
A-3. 

5. Jet Propulsion Laboratory emphcracris tape. This tape 
is accessed to obtain the positions of moon and sun. 

Radio Star Locations. There arc five radio stars at 
both 156 MHz and 400 Mllz (in Section 5.3). 

Tljc right, ascension, declination and noise power for 
cacl) radio star is sliown i?i Section A- 4, 
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A-1 RAE-B TRACKING STATIONS 
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■ APPHKmX P. 

l)P.TATU:i) COMi’HTA'l lOXAL PUOCpmiUliS 
AND I'l’.nGPAM Pi.OV.- CHARTS 


Tbir. 
s ir. p 3 i f i o il 
di i! v.lii 

brariches . 


f.cc.t^on tlcscv jbcs each processing step in the 
prc-:r.ra and is fol lov.td by an equation flow 
ell shov.’s ll'.o cor.putationnl Tlo\.' and logic 


Find o: (right ascension) and 6 (doclinaticn) of moon 
and sun at inertial time from the 
JPI/ cphc!.ieris tape, where the (x,y,z) of Sun and 
(x,y,z) of Moon (coordinator referenced to equinox; 

arc given. 

Then 


a = tan — 
X 


9 



\ 


and 


module^ (o: , 2u) 



l);;o x,y ,2 coordiii.-it c- 
the;', ca 1 c.iil to Llic a 


to co'.n'.^uto the unit vector, 
Ic bctv.’cc'U Moon and Sun 


I'or:' cart li roon unit vector and earth-sun unit 
V i- c t o r 1 . 

b 


moon-earth" 3 . 


m 




laj 


^mk 


in ' 


. 1 


m 


9 •> 


»OC^+Y 

m m 


m 


m. 

3 


A ^+Y ^^+2^ 
m m 


m 


"'l: 


7 x 


_ 2 

( '-+Y +2 

m m m 


Note : 
• ^ 


Assumes antenna borcsight is at moon's center. 


sun- earth 


SI 


hj 


r»l: 


jy ^+y 




z 


5 


, J ^ 


b. Calculate t)>e anpjc botu'ccii Moon and Sun. Use 

two unit vectors T and T comiuitc cross product 
• ,s m ' ' 

Sin'^ IT,. I) 

S. Vi.slbility Test 

a. For eadi station, use the Crocnv;ich Sidereal Tiric to 

get the position in the time of observation. 


Let X be the east longitude, & be the right ascension 
of Greenwich at the tirae 

Julian Date JD = 2433283.423+565.2422(1973-1950} 

. + Day of the year interested, 

„ . ^ JD-2415020.0 

Fraction of Julian Century T = 

36525 


®gQ 99*’,6909833+56000® 7689T+0“.0003S70ET^ 

®p "" 0 . 2.S06S44 7r(Min.} 

a (right asccn.sion) =0 + X 

E 

a = AMOD (a, 6.2S3185) 


6 (declination) . = Tati tuclc 


00 


I 



For cocli still i Oil, the unit vccLoi is 


cos p cos 6 
sin a cos 6 
si n 6 

Then calculates the Dot Product of two unit vector 

Tt , • Tms • T 

T’ m -^T in 

If the result of dot product is negative, set temper- 
ature equal to zero and continue to tost next station. 
Otherwise, select antenna corresponding to station and 
the temperature of the antenna for that tine point is 
t cjii)H:raturc of: station. 

Sun t cnii'» (.vr a t u V c ^ ) oal cu 1 ci t i on 

If sin"^ Igl) 1 - lobe beam width, perform 

coii^]nit at ion. 


Hi 


11 


1, 


Tk 


^Ti 

Hj 


'Tk 


r>] 



- 0 • - angular dir.;:'.;: ter of sun's apparent 

temperature model {- • 66®). 

= full bcar.iv.i dth of syminctri cal antcnr.a 
jiiain Icbc. • 

= quiet sun ideal model. 
s;6 X 10^ K® for 400 MHz: 

8 X 10^ K° for 155 MHz. 


If any antenna first side lobe v.ould contain sun 
Y » 1/.2 0^ (PRIMARY) + 6^^ (SECONDARY) 

0^ (PRIMy\RY) == full bermv.’idth of main lobe. 

0^ (SECONDARY) - full bcamv.’idth of first side lobe. 



Perform computation as before except substituting 
0^ (S]-CONJ)ARY) ->■ (PRIMARY) and multiply by 



so, the; 1 Ion is 


2 

Ts ” Tb 


pov.cr R;:in 
10 


Pov.cr gain (HB) : 

(Avera'^c first 
side loiic gain) 


-IS. 5 -- in VHP lUUK system , 136 Milt 
-16. 25-- in 85 FT antenna, 4pOM!lz 
-30.0 -- in 40 FT antenna, 400>5Hz 


Otherwise , 


set equals to zero. 


.S. For cacli Radio Sta-r coir-putc the position .at prediction 

time and if the angle between the Moon and Radin 
Star is greater than antenna half boamwidth, set 
the temperature du.c to Radio Star cciual to 

zero, otlicr'A'isc , compute the temperature . 


llLilizing tlie right ascension and declination, compute 
the unit vector. Computation’ is the same as in 3.b. 





C;i] dilate the cre?5 product of tlic two unit vectors; if 

^ I f ^ — 7 iQbc- beam width, cor.ijiutc tlic 

temperature due to the Radio Star. Otherwise, set tlic 
tciiiperaturc ccjual to zero, then test next Radio Star. 

Core/iutation of Radio Star temperature 

1 *^r 

"^STAR ~ 7 ^ (Kelvin decrees') 

2 . 47T K i> . 

is the receiving antenna gain above isotropic. (Dii) 

X is the wave lengtli (.meter), 

observed radio flux density (VAatt m'-^Kz"'^) 

K Boltcmann's constant.' (1 .38x10’^^ J/K) 

After test all five radio starts, acid each temperature 
to get tlie tot;il tenpera ture . 


Compute Brigiitness Temperature 

From reference 2, the follov.’ing formulation is given; 



T (e,<») G (0,<)) sin 0809<j) 
G (0,^0 sin G905(f. 


G = elevation betv.’ecn boresi[’,ht axis and the 
o 

horizon (deg.) 

'l’(0>'i') is noise temperature distribution of the 

galax)' (excluding the sun and predominant radio stais) 
obtained from input data. 

G(0,<0 is the lossless antenna gain distribution. 

0 and ({• are the orientation angles defining the 
positioii of a radial surface element within the 
celestial hemisphere. 

In this program, the numerical method used to solve this 
equation is as follows: 

From Moon’s right ascension (a) and declination (6) 
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